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Abstract 

In this letter we explore the dependence on the gauge fixing condition of several quantities in the 
U{\) Higgs model at finite temperature and chemical potential. We compute the effective potential 
at the one loop level, using a gauge fixing condition that depends on \x and which allows to decouple 
the contributions of the different fields in the model. In this way we get the mass spectrum and 
the characterization of the phase transition, pointing out in each case how these quantities depend 
on the gauge fixing parameter £. When /i vanishes, we agree with previous results if £ = 0. The 
gauge dependence problem is also analyzed from the perspective of the Nielsen identities. 
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I. INTRODUCTION 



In spite of being a very simple case, the U(l) Higgs model acquires a quite complicated 
structure when we go to the scenario of finite temperature and chemical potential fi. For 
example, since the paper by Bernard [1], it is well known that at finite temperature the 
ghost fields cannot be factorized in a trivial way from the functional integral in the U(l) 
case, as it does occurs at zero temperature. 

The gauge fixingcondition at finite temperature is a long standing problem in field theory. 
Dolan and Jackiw 9j explored the gauge dependence of several quantities, in particular the 
critical temperature from the effective potential. Keeping in mind the validity of the one 
loop calculation up to order 1//3 2 , it is possible to restore the gauge invariance by expanding 
the effective potential up to this order. 

The relationship between the gauge fixing condition and the temperature for the phase 
transition in the electroweak theory has been explored in 10J. The analysis of the phase 
transition at finite T, by the Landau method has been discussed in llj. In particular, in 
the last article, the authors found a critical temperature that does not depend on the gauge 
parameter. However, this is due to the fact that a gauge condition has been chosen, which 
couples the Higgs with the gauge field. As we will see in this article, if we decouple the 
Higgs from the gauge field, which is another possibility for fixing the gauge condition, the 
critical temperature still depends on the gauge parameter £. 

In the literature, people usually chooses the £ = 1 Feynman gauge, or £ = Landau 
gauge, claiming convenience reasons. When discussing the effective potential in the frame of 
the Weinberg-Salam model, the unitary gauge is excluded because it leads to an erroneous 
result for the critical temperature and the pressure of the system [2']. It seems that as soon 
as temperature is turned on. the gauge invariance of the theory becomes a subtle subject. 
In 3)], the authors emphasize that only periodic gauge functions (A p ) are allowed, i.e. those 
that satisfy A p (x M + if3u^) = A p (x M ). In fact, the gauge fixing procedure and the gauge 
invariance of theories at finite temperature and density has not yet been fully clarified. In 
particular, if we expand the effective potential up to the order l/(3 2 following jsj], the gauge 
dependence is not removed when finite chemical potential is present. 

In this letter, we will consider the U(l) Higgs model in the presence of finite temperature 
and chemical potential. First, we calculate the effective potential at the one loop level, using 
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a gauge fixing condition that depends explicitly on the chemical potential, and which was 
used previously in a different context 7|,|8j]. In this way, we get non trivial dispersion relations 
for the different particles of our model. The role of the ghost fields is analyzed in detail for 
different cases. We compare our results, when /j, — 0, with previous calculations [3], been in 
agreement with the results by Kapusta when £ — > 0. We calculate the critical temperature 
T c associated to the restoration of the U(l) symmetry, emphasizing its gauge dependence. 
The relationship between the chemical potential and the gauge invariance problem is also 
analyzed. 



II. THE 17(1) HIGGS MODEL 

Our model is described by the following Lagrangian 

C = (D^YD^ - m 2 (0*0) - A(0*0) 2 - \F^F^ (1) 

where 



D p <j> = (d p + ieA p )<P, (2) 
F pv = d p A v - hAp. (3) 

The set {<ft, <j)*}, represents a pair of complex scalar fields. We will start in the broken phase 
(m 2 < 0), where we have a real non- vanishing vacuum expectation value u for the <fi field. 
The local gauge transformations U that leave our Lagrangian invariant are the following 



(A^(x)) u = A"(x) + d£A(x), (4) 

4> U = * u -v, 

= 4>{x) — ieA(x')4>(x') — v. 

The partition function is given by 

Z = N{f3) I D(j)D(f)*T[DA p exp I dr [ d 3 xC eff , (5) 
./Per. Jo J 
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where Per in the functional integral indicates that we have to integrate over periodic eu- 
clidean fields configuration in the interval (0,/3). 

The Lagrangian £ e ff i n the previous equation, incorporates the chemical potential 
through the recipe of considering fi as the zero component of a constant external gauge 

field y, 



C 



eff 



■(d p - ieA p + ifi5 p0 )4>*(dp + ieA p - i/i5 pO )0 
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(6) 



III. ONE LOOP EFFECTIVE POTENTIAL 



In order to calculate the effective potential, it is more convenient to express the complex 
6 field linear combination of real fields 4>i an d (p2, such that (f> = ^(</>i + i<j>2)- 4>i 
acquires the non vanishing expectation value v. 

We will expand the action up to second order in powers of the fields, using a Lagrangian 
C j that incorporates external sources, J\, J 2 and J M , 



J 1 , h 



Jn 



In this way we get 



S= j P dr I d 3 x£ T = + «S« + S (2) + 



where 



S (o) 



£ 2 
2 



m 2 n A a 

— v v 

2 4 



(3VL + Jiz/, 



and 



(7) 



(8) 



(9) 



= J dr J d 3 x.{—ifiv5podp(j)2 + ^ 2 ^4>i — eA p 5 p0 ifi 



(10) 



v 2 — m 2 u(j) 1 — Xu 3 (f)i 

J\ , J2 , Jp 

+ ^ + + w 

By shifting 4>i — > 4>i + v i we can suppress terms proportional to 4>i in since the 

quantum fields must have a zero vacuum expectation value. Choosing 
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i = z/(m 2 -fi 2 + Xu 2 ), 



and 



(11) 



Ja ■ 2 

we can eliminate the terms proportional to <pi and A in 
The <fv > can be written quadratic form 



(12) 



( ^ \ 



\ Aa J 



where M is a matrix given by 

' -dpd p - i? + m 2 + 3\v 2 



2i[id p 5 



\ 



2i^5 0a v 



dpdp — [i l + m 2 + \v 



-evd a 



2ifidopU 
evdp 
e 2 5 pa u 2 + B pa j 



The operator B pa is defined as 



(13) 



(14) 



B pa = -d x d x 5 pa + d p d a . (15) 

In the previous matrix we have undesirable mixing terms which couples (j) with A. These 
terms can be eliminated through an appropriate gauge fixing condition given by 

F = (dp + i2^ pQ )Ap - iefrifa + ifa). (16) 

This gauge fixing condition was introduced previously [7] to calculate the contribution of 
each field to the the effective potential of the Weinberg-Salam model, and as a background 
field that maintains explicitly the gauge invariance 

From this condition, the Faddeev-Popov Lagrangian, at finite fi can be read as 



Cfp = d p fj{d p + 2tfi5 p0 ) V - ^e 2 u 2 V (17) 

Finally, our total lagrangian is written as a sum of the different contributions Cr = 
Cj + Cqf + £fp- The ghost Lagrangian plays a important role when counting the degrees 



of freedom of the model, but since they do not couple to the physical fields, their role will 
be discussed in the next section. The action term is defined by a new matrix M, 



/ 



V 



d p d p - /i 2 



\ 







2ijj,d p 5po 
-dpdp — fi 2 + ml 



e 2 5 pa is 2 + B pa j 



(18) 



where we have introduced effective gauge dependent masses mi, m 2 as 



mf = m 2 + 3Xu 2 + eV£, (19) 
m \ = m 2 + \v 2 + eV£. (20) 

Note that both mi and m 2 depend on the gauge parameter £ this means, for an arbitrary 
value of £ there are no Goldstone boson in this model. We will come back to this point later. 
The operator B pa is the extension of (fl~5l) including gauge dependent terms 

B pa = ~d X d X 5pa + dpd a (l ~ + V^W'- ( 21 ) 

Finally, the effective potential is given by 
with 



Z 1 = DfaDfoDAe^ 2 ', (23) 
./Per. 

where Q is the spatial volume. 

In the momentum space, In Z\ can be written as a sum over the Matsubara frequencies 

uj n = 27in/j3 

InZi = -§ En I $js In {det pa ([£; 2 + eV]^ - fc p fc a (l " T 1 ) - V^pokor 1 )} 
-§ E„ / ^ ^ {k* + k 2 {m 2 + m 2 - 2^ 2 ) + 4^ 2 n + (m 2 - ^{m 2 - /, 2 )} , 



(24) 
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with k 2 = u 2 + k 2 and where det pQ refers to the determinant in the Lorentz indices. The 
polynomials in the equation above, since they are quadratic in u 2 , can be factorized such 
that 

In^l = -§ En I j£$r ln{(«* + + 
-§ E„ / iftr M( W 2 + y 2 )(^ 2 + y 2 2 )} + ln{(^ 2 + y 2 )( W 2 + y 2 )}, 



(25) 



Using the identity 



£ ln(^ 2 +X 2 )=pX + 2 ln(l - e"^), (26) 

n 

we can decompose the effective potential as a sum of terms associated to the different fields 
in the model. In other words we have diagonalized the effective potential. Then 

V eff = V* + V A , (27) 

where V^a is given by 

V --]J^— + P ' (28) 

which is valid for each component of the fields. 

In this way we can identify and yi (i = l.A) with the energy spectra of these pseudo- 
particles 



x 1 = m-i 



I + k 2 - Az/ 2 + /i 2 + ./4/i 2 (m 2 + k 2 - Az/ 2 ) + AV 4 , 



x\ = m 2 + k 2 + Az/ 2 + fi 2 - \j^ 2 {m\ + k 2 + Az/ 2 ) + A 2 zA (29) 
Notice that if v = and // ^ 0, (129!) simplifies to the well known expression 



Xl >2 = El = Vk 2 + m 2 ± /i. (30) 

Here, mi = m 2 = m and there is no dependence on £. If /x = and z/ 7^ 0, xi and £2 depend 
on £. 

For the gauge fields we find two non-anomalous massive "photonic" states which obeys 
the following dispersion relation 



yi = y 2 = E J = Vk 2 + eV 2 . (31) 

Notice that x\ and x<i describe the possible energy values for 4>i and fa. 
The pathological energies associated to the other two "photonic" degrees of freedom (1251) 
are given, in terms of e = 1 — £, by 



2/ 2 =k 2 + ^(1 + + 2^ 2 - iy^W+W^lWW). (32) 

We call these energies pathologicals, since they do depend explicitly on the gauge param- 
eter £. 

Repeating the same previous analysis, if v — and /i / 0, with £ = 1, we recover a 
photon with energy y 4 = |&|. 

On the other side, in the broken phase when \i — 0, we have 1/3 = Vk 2 + e 2 i/ 2 an d 
?/ 4 = A/k 2 + eV 2 £, i. e . we have three massive degrees of freedom (3/1,3/2 an d 1/3) for the 
photon and a gauge dependent contribution (j/ 4 ). Now taking /i 7^ and £ = 1, we recover 
the three degrees of freedom for the massive photon: yi, 1/2 and 3/4. 

As we will see in the next section, the Faddeev-Popov ghost fields will compensate these 
anomalous dispersion relations. In particular for £ = 1, the situation is completely clear. 
This probably is one of the reasons why people prefer the Feynman Gauge £ = 1. 



IV. GHOSTS AND EFFECTIVE POTENTIAL 



In order to get a proper counting of the degrees of freedom, we have to introduce the 
Faddeev-Popov ghost fields. This point is completely different with respect to what occurs in 
the U(l) theory at zero temperature, where the ghost can be factorized and do not contribute 
to the effective potential. 

The functional integration over the ghost fields, represented by the complex Grassmann 
variables (r),fj), leaves 



DrjDrj exp 




d 3 xr] - 



5_F_ 
IX 



n 



ln(u;„ + k 2 + eV£ + 2ifjuj. 



(33) 



If we write the argument of the logarithm in the polar representation 



pe* = u 2 n + k 2 + eVf + 2ipu n , (34) 

with 



P 



vVn + k 2 + e 2 z/ 2 ) 2 + 4/A4, 
4u n p 



arctan 



uo 2 + k 2 + e 2 v 2 J ' 

we see that the sum over vanishes since (p(uj n ) = —ip(—uj n ). Factorizing p 2 
4)^1 + 4), we find 



(35) 
(36) 



Zl,2 



E v ± = Vk 2 + e 2 z/ 2 £ + p 2 ± p. 



(37) 



If p — 0, it is interesting to notice that one of the ghost fields cancels the contribution 
from the nonphysical photon y% 

Using EBl we can decompose the effective potential such that V e ff = V tre e + V^ + Va + V v , 
where the different contributions are 

• Tree level contribution 



V t 



tree 



m 



A 



if + -v 
2 4 



contribution 



(38) 



f d 3 k (xi(Jfe) +x 2 (k)) ln(l- e -^ 1 ( fe ))(l-e-^ 2 ( fe )) 

Massive photons contribution 

d 3 k (yi(k) + y 2 {k)) ln(l - e^ l(fc) )(l - e -^ 2(A;) ) 



(27T) 



+ 



/3 



+ Vk 2 + e 2 u 2 + - ln(l - e -W+<?u^ 
P 



Ghost contribution 



d 3 k (zi(k) + z 2 (k)) e~^ k) ){l- e-^ k) ) 

Vn — / ~ ~ ~ 



(27T) 







(39) 



(40) 



(41) 
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V. HIGH TEMPERATURE EXPANSION 



Our previous results for the effective potential are interesting since we decouple the 
contributions for the different fields. However, the integrals cannot be calculated analytically, 
and therefore it is appealing to carry on a high temperature expansion of the effective 
potential, which will allow us to compare our expressions with well known results from the 
literature. 

For the one loop effective potential in a high temperature expansion we find 



m 



A 



— v - —V* + -Z/* 

2 2 4 



Veff 

The vacuum is defined by 



2- 2 7 1 / -' , , , ., 

+ — {m 2 + 2\v 2 

45 12 ^ 



(42) 



(fi 2 + \m\ 2 ) 



1 _ 11 

1- - T 2 



para T <T C 



MT)Ln={ L icJ ~ (43) 

0, para T > T c . 

= 12(fi 2 + \m\ 2 ) 

c 4A + 3e 2 + £e 2 ' 1 ' 

where T c is the critical temperature where the symmetry is restored. 

The critical reader at this moment could be surprised that the critical temperature, which 
is in principle an observable quantity, depends on the gauge parameter. Notice, however, 
that for finite chemical potential it is well known that the number of Goldstone bosons is 
lesser than the usual prediction ( 12]). This is in agreement with our results. In our case, 
both higgs fields acquire a mass. On the other hand, if fi vanishes, we know that we must 
recover one Goldstone boson. This is not possible, unless that £ also vanishes. In this case 

ri r 

the i?£ gauge becomes the Lorentz gauge. Our results is then in agreement with [2j and [10| . 
when [i = 0, in the gauge defined as £ = (Landau gauge). So this gauge seems to be the 
more appropriate to extend the calculations to the scenario with finite mu. 

For finite /x, it seems that the one loop calculation of the effective potential is not enough 
for having a gauge independent result for the critical temperature. 

In general, from previous work [13J] we know that the determination of critical values of 
parameters associated to phase transitions requires to go beyond the one loop approximation 
through an appropriate re-summation. 
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VI. SYMMETRY RESTORATION, PHASE TRANSITION AND NIELSEN IDEN- 
TITIES 

As it is well known, the occurrence of the phase transition can be inferred from the be- 
havior of the isothermals of the effective potential. The pressure is defined as P = P-§y In 
which implies P = —V e ff, since the effective potential corresponds to the thermodynamical 
potential of the system. 

We found for the pressure. 



P < = — + 12 + 4A 1 1 "^) ' (45) 



2?r 2 T 4 \m\ 2 - u 2 f- 1 
> 45 12 y ' 

where P < (P>) corresponds to the broken (symmetric) phase. Here we have the same 
difficulty, a gauge dependence, we found previously for the critical temperature. 

The pressure and the entropy are continuous at T c . However the specific heat has a 
discontinuity which depends on £ only through T c . This means that this is a second order 
phase transition. In fact, the behavior of the specific heat confirms this picture. 

The gauge dependence problem can also be analyzed from the perspective of the Nielsen 



Identities. This method [I4j], related to the BRST symmetry transformations, is a procedure 
that allows to search for possible gauge dependence of physical quantities which are related 
to an explicitly gauge dependent effective action. If these identities are satisfied, we may 
have confidence that the results for the physical magnitudes will be gauge independent. The 
Nielsen identities are still valid for finite temperature [jjj]. An interesting discussion of the 
Nielsen identities for the generalized gauge in the abelian Higgs model can be found in 
16l | . We will follow the strategy of this article to explore the validity of the Nielsen identities 
for our case. 



3- 

( 

The Nielsen identities arise from the following identity 



F_dF_ 



(47) 



where we use the notation 
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(O) = e r J DfaDrjDfjO exp (-S F + J d^x^-fa - Pi)) , (48) 

and where i], fj are the ghost fields, F is the gauge fixing condition, which in our case is 
given by eq. (TIE]) . T is the effective action, £ is the gauge parameter and Aj is the BRST 
transformation of the scalar fields <\>\ and 02 which is given by S(pi= Ait] with Aj = (—02, 0i)- 

In the finite temperature scenario, we have to keep in mind that the integrals in the 
previous equation, when going to the momentum space must be handle (d=4) as 

/ «** -E/(05. <«) 
where we sum over Matsubara frecuencies u n = 2im/i3. 

Assuming for the classical components of the fields 0i, 2 , 0i = v and 2 = 0, the Nielsen 
identities can be expressed as 



where C is given by 



c = \ fd d x f d d y (MyMy)v(x) [(d p + z2/^ )A> + ie£u(<l>i + # 2 ).]) , (51) 

n 

and must be evaluated in a perturbative expansion, see [16]. 

We explored the validity of the Nielsen identities for the following cases 

• [A — 0, v — 0: It is very easy to see that the equation (!50i) is valid for any value of £, 
in ths case, confirming [ijj] the validity of the Nielsen identities for finite temperature. 

• fi = 0, v In this case the Nielsen identities are valid only if we take £ = and if 
we evaluate the derivatives in v = v m i ni i.e. in the classical value. Note that £ = is 
demanded by the Goldstone theorem, which implies rri2 = 

• jj, 0, v = and fi ^ 0, v ^ 0. In these cases the Nielsen identities are not valid 
for any possible value of £. In fact, the only non vanishing contribution to the left 
hand side of the Nielsen identities, according to eq. ( 1501) is given by the "photonic" 
dispersion relations, eq. (1321) . for the fi = 0, v ^ case. When /i ^ 0, v ^ 0, all 
degrees of freedom contribute to the derivatives that appear on the Nielsen identities. 
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However, in both cases the diverge when £ goes to zero. This means that we are 
not allowed to take the £ = g auge, where the Goldstone theorem is realized, in 



agreement with the result by 12]. 



To summarize we can see that the \x dependent gauge fixing condition we used here, and 
which was valuable for diagonalizing the effective potential in the Weinberg-Salam model 
for £ = 1, p|, turn out to be unnatural when looking for gauge invariant predictions for 
physical quantities in the abelian Higgs model. 

The gauge dependence problem is related to two different aspects: 

• The necessity of having gauge transformation which are periodic in the temporal di- 
rection, with period j3 and 

• the occurrence of finite \i. 

As a conclusion, we would like to remark that the problem of the gauge fixing at finite 
temperature and/or density has not yet been solved. In the U(l) Higgs model we have 
explored in detail how this gauge dependence propagates through several physical relevant 
quantities. Nevertheless, when \i vanishes, we agree with previous results in the literature. 

It can be shown that /i can be incorporated as a boundary condition for the Green function 
of fields confined to a finite spatial region, as it occurs in chiral bag models ^, [f| . Eventually, 
the construction of the theory should start by taking the theory in a bounded space region, 
with adequate //-dependent boundary condition. When taking then the limit where we go 
into the whole space, probably a remanent of the non-trivial boundary conditions will remain 
as a non trivial dependence on the gauge fixing condition. This point will be explored in a 
future work. 
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